arXiv: 1505.07669vl [astro-ph.GA] 28 May 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 29th May 2015 (MN KTeX style file v2.2) 


Imaging the cold molecular gas in SDSS J1148 + 5251 at z = 6.4 

Irina 1. Stefan/* * Chris L. Carilli/’^ Jeff Wagg/’^ Fabian Walter/ Dominik A. Riechers/ 
Frank Bertoldi/ David A. Green/ Xiaohui Fan/ Karl Menten/ Ran Wang^’^ 

^ Cavendish Laboratory, 19 J. J. Thomson Ave., Cambridge, CB3 OHE, UK 
^National Radio Astronomy Observatory, Socorro, NM, USA 

^Square Kilometre Array Organization, Jodrell Bank Observatory, Lower Withington, Macclesfield, Cheshire, SKll 9DL, UK 

‘^Max-Planck Institut fiir Astronomie, Kbnigstuhl 17, D-69117 Heidelberg, Germany 

^Department of Astronomy, Cornell University, 220 Space Sciences Building, Ithaca, NY 14853, USA 

^Argelander-lnstitut fur Astronomic, Universitdt Bonn, Auf dem Hiigel 71, Bonn, D-53121, Germany 

''steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 

^ Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

^Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China 


29th May 2015 


ABSTRACT 

We present Karl G. Jansky Very Large Array (VLA) observations of the CO (J = 2 1) 

line emission towards the z = 6.419 quasar SDSS J114816.64 + 525150.3 (J1148 + 5251). 
The molecular gas is found to be marginally resolved with a major axis of Oi'O (consistent 
with previous size measurements of the CO (J - 7 ^ 6) emission). We observe tentative 
evidence for extended line emission towards the south west on a scale of ~ T.'4, but this is 
only detected at 3.3cr signihcance and should be conhrmed. The position of the molecular 
emission region is in excellent agreement with previous detections of low frequency radio 
continuum emission as well as [C n] line and thermal dust continuum emission. These CO 
(J = 2 1) observations provide an anchor for the low excitation part of the molecular line 

SED. We find no evidence for extended low excitation component, neither in the spectral line 
energy distribution nor the image. We fit a single kinetic gas temperature model of 50 K. We 
revisit the gas and dynamical masses in light of this new detection of a low order transition of 
CO, and confirm previous findings that there is no extended reservoir of cold molecular gas 
in J1148 + 5251, and that the source departs substantially from the low z relationship between 
black hole mass and bulge mass. Hence, the characteristics of J1148 + 5251 at z = 6.419 
are very similar to z ~ 2 quasars, in the lack of a diffuse cold gas reservoir and kpc-size 
compactness of the star forming region. 

Key words: galaxies: active - galaxies: formation - galaxies: high-redshift - galaxies: indi¬ 
vidual: SDSS J114816.64 + 525150.3 - cosmology: observations - radio lines: galaxies 


1 INTRODUCTION 

Studies of molecular lines towards high redshift galaxies provide 
powerful diagnostics on the properties of the dense interstellar me¬ 
dium (ISM) in the early Universe. As cool gas is the fuel for star 
formation, such studies are probes of galaxy formation and evolu¬ 
tion and help us to understand the origins of today’s galaxies (for a 
recent review, see |Carilli & Walter|2013| l. 

Atz = 6.419, SDSS J114816.64-t525150.3 (hereafter J1148 + 
5251) is the most distant quasar discovered in the Sloan Digital 
Sky Survey (SDSS, [Fan et al.|2003[ l, and remains among the most 
distant known. It is a radio quiet quasar ^Carilli et al.|2004| jPet-j 
|ric et al.]|2003^ with an estimated supermassive black hole mass 

of 3 X 10^ Mq accreting at the Eddington limit, as found from 
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Mg II emission lines jWillott, McLure & Jarvis|2003| ). Optical spec¬ 
tra of the source show a saturated Gunn-Peterson trough due to a 
partially neutral intergalactic medium, providing some of the first 
evidence from the earliest epochs of galaxy formation and cos¬ 
mic reionization at z > 6 ( [White et al.||2005t . J1148 + 5251 was 
the first z > 6 quasar to be detected as a hyperluminous infrared 
galaxy (HyLIRG) (infrared luminosity > 10'^ Lq, 

|2003a^ through millimetre bolometer observations. It also remains 
the most distant detection of molecular gas to date, with multiple 
detections of CO rotational transition lines from 7 = 3 —> 2 up 


to 7 = 7 —> 6 ( [Walter et al.|2003||Bertoldi et al.|2003b[|Riechers| 


[et al.|2009 
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leus (AGN) and extended on scales of 5 kpc ( [Walter et al.| 2004l 
[Riechers et al.|200^ . 

[Cii] observations trace a bright 1.5 kpc atomic gas region 
with active star formation apparently taking place at an estimated 
star-formation rate (SFR) surface density of ~ 1000 Mq yr“' kpc“^ 
( [Walter et al.|2009| (. Based on [C n], [Cicone et al.[ ( [2015[ l suggest the 
presence of a high velocity complex outflow extending to projected 
radii of ~ 30 kpc. J1148 -(- 5251 has a far-infrared (FIR) luminosity 
of 2.7 X lO'^ Lq Leipski et al.|2013} , co rresponding to ther mal 
emission from 4.2 x 10“ Mq of warm dust jBeelen et al.|2Q06^ . 

J1148-(-5251 has served as the archetype for extreme starburst- 
AGN at high redshifts, representing the early, coeval formation of 
supermassive black holes and their stellar host galaxies during the 
epoch of reionization. 

In this paper, we present observations with the Karl G. Jansky 
Very Large Array (VLA) of the CO (/ = 2 —> 1) line emission 
in J1148-(-5251. These observations have a number of significant 
advantages over the previous VLA observations. Low order CO 
transitions are critical to constrain the total molecular gas mass in 
galaxies, and in particular, to look for massive reservoirs of cold gas 
fuelling star formation ([Riechers et al.|2006||2011a|b||Ivison et al.| 
[2011[ |. The lowest previously detected transition from J1148-(-5251 
was CO (y = 3 —> 2) [Walter et al.[[2003) , with rather shal¬ 
low limits on the CO (7 = 1 —> 0) line emission ( [Bertoldi et al.[ 
[2003b[ (. CO (y = 3 —» 2) can be substantially sub-thermally ex¬ 
cited, even in starburst galaxies, with a typical brightness temper¬ 
ature ratio of 0.66 for submillimetre galaxies (SMGs), and 0.6 for 
high z main sequence galaxies ( [Carilli & Walter[2013[ . Conversely, 
the CO (y = 2 —> 1) to CO (y = 1 —» 0) ratio is almost always 
thermal (i.e. the brightness temperature is constant due to thermal- 
ized emission), even in high-z main sequence galaxies, with ratios 
of 0.85 in SMGs and 0.9 even in high z main sequence galaxies. 
Hence, CO (y = 2 —> 1) provides an excellent route to study cool 
gas reservoirs, being four times stronger (in erg s“* Hz“*) than CO 
(y = 1 —> 0), but still tracing the colder gas ( [Daddi et al.[2014[ l. 

We use the capabilities of the VLA to image the full emission 
line from 11148 + 5251. The previous CO (J = 3 2) observa¬ 

tions used the old VLA correlator, which had limited bandwidth 
and channelization, leading to an observation that missed the line 
wings beyond ±150 km s“*, and had only 12 channels across the 
centre of the line. 

And lastly, we obtain sensitive observations of the continuum 
emission at a rest frame frequency of 230 GHz, corresponding to 
an emission regime that would be dominated by the cold dust. 

We adopt a standard ACDM cosmological model with Hq = 


67.3 km s ' Mpc ', = 0.685, Om = 0.315 (Planck Collabora- 


[tion et al.|2014[ l. At a redshift of z = 6.419, this implies an age of 
the universe of 851 Myr and a scale of 5.623 kpc arcsec^*. 


2 OBSERVATIONS 

We use the VLA in C and D configurations to observe the CO (J = 
2 —> 1) (vrest = 230.538 GHz) emission line towards J1148±525L 
The line is redshifted to v = 31.074 GHz at z = 6.419. Observations 
in the D configuration were performed between November 2011 
and January 2012 for a total of ~ 15 hours on source. Observations 
in the C configuration followed between January and April 2012 
for a total of ~ 9 hours on source. Since these observations were 
taken as JVLA early science during the initial commissioning of 
the new correlator, we chose to observe topocentric frequencies. 
Furthermore, the observations span a six month range that involved 


Table 1. Noise and resolution of a number of maps obtained with the VLA 
data presented in this paper. 


Arrays 

used 

‘robust’ 

parameter 

rms per channel 
(19.309 km s-') 

Restoring 

beam 

c 

1 

95 i-iJy beam“* 

0'.'76 X O'.'69 

D 

1 

68 i-iJy beam“* 

2'.'53 X 2'.'28 

C&D 

-1 

126 [iJy beam“^ 

a'55 X a'52 

C&D 

0 

73 i-Jy beam“* 

a'82 X 0'.'75 

C&D 

1 

56 i-iJy beam“* 

1"75 X T.'62 


major changes to the correlator between C and D configurations. 
Hence, we did not attempt to doppler track, but have summed the 
topocentric spectra. This leads to at most a 20 km s“* smearing of 
spectral features, which is comparable to a single channel in our 
final analysis. This smearing will have no effect on our subsequent 
analysis. 

The standard calibrator 3C286 was used as a primary flux and 
bandpass calibrator and the nearby source J1153±4931 was ob¬ 
served approximately every five and a half minutes for pointing, 
secondary amplitude and phase calibration. 

These observations were taken during the commissioning 
phase of the VLA correlator, and hence correlator capabilities were 
evolving rapidly. The D array observations were obtained when full 
channelization (64 channels) was available, but only for a 128 MHz 
bandwidth. This bandwidth is easily adequate to cover any reason¬ 
able line emission (±600 km s“*), but does not allow for a sensitive 
continuum measurement. The later C array observations had the 
same line coverage, but also had close to 1 GHz of bandwidth off¬ 
line, providing a sensitive continuum observation. 

The data were calibrated, imaged and analysed using standard 
techniques in NRAO’s CASA packag^ ( [McMullin et al.|[2007 1 . 
Some difficulties were encountered due to the presence of two 
bright sources near the half power position of the primary beam. 
Small inherent pointing inaccuracies resulted in these sources ap¬ 
pearing as time-variant leading to significant residual sidelobes 
even after deconvolution. To solve this issue we performed con¬ 
tinuum subtraction in the Mv-plane using the outer most channels (4 
to 9 and 57 to 62) of the 128 MHz band centred on the emission 
line (channels 23 to 43). 

We imaged the C-array data separately using a close-to- 
‘naturaT weighting scheme characterized by a ‘robust’ parameter 
( [Briggs|1995) of 1. We then repeated the process for the D-array 
data only. We also created a combined data set from the C and 
D-array data and imaged it with different weightings to explore 
spatial resolution versus signal-to-noise optimization. A close-to- 
‘uniform’ weighting scheme (robust = -1 in CASA) was used to 
achieve higher spatial resolution, a close-to-‘naturaT robust = 1 
was used to obtain higher sensitivity, and a good trade-off between 
resolution and sensitivity was achieved with a middle-ground ro¬ 
bust = 0. This way we produced maps with a large range in resolu¬ 
tion, as shown in Table [T] 


3 RESULTS 

The field covered by these observations also includes two bright 
radio sources [Becker, White & Helfand|1995[ l, one ~ 35" to the 

* The Common Astronomy Software Applications package, 
http://casa.nrao.edu 
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J2000 Right Ascension 

Figure 1. Continuum map of the C-aiTay data. The rms is 6.4 pJy. Two 
lower redshift sources to the north west and the south east of J1148 + 5251 
dominate the emission all over the map. The cubehelix colour scheme 
jGreen|201ll was used. 



Figure 2. Spectrum of the CO (7 = 2 —» 1) line emission towards J1148 + 
5251 obtained with VLA D array observations at a resolution of 19.335 
km s“' (2 MHz). The blue line shows a Gaussian fit to the data. The velocity 
scale is relative to the redshifted CO (7 = 2 —> 1) frequency at z = 6.419. 


NE of J1148 + 5251 and one ~ 40" to the south west, both at 
much lower redshift ([ Becker, White & Helfand|1995^ . They have 
been mentioned in previous studies of J1148+5251 (e.g. [Bertoldij 
jet al.|2003al l. Figure[^is a continuum map of the C-array data only, 
showing these sources. We measure a non-primary beam corrected 
flux density of 2.79 ± 0.05 mjy for the north eastern source and of 
2.98 ± 0.06 mJy for the south western. We perform uv-continuum 
subtraction to remove the effects of these sources from subsequent 
analysis. 

We clearly detect the CO (7 = 2 —> 1) emission towards 
J1148 + 5251. Figure]^ shows the CO (7 = 2 —> 1) emission 
spectrum obtained from the D-array data. Fitting a Gaussian to the 
emission gives a peak intensity of ly = 0.298 ± 0.022 mJy beam^' 
and a FWHM of 298 ± 26 km s“' centred on -5 ± 15 km s“*. This 
implies an integrated line flux of 94.6 ± 7.7 mJy bearn^* km s“' and 
a line luminosity of = 3.2 ± 0.3 x 10'° K km s“' pc^. 

Constraining the source size in CO (7 = 2 —> 1) emission is 



Figure 3. Spectra of the CO (7 = 2^1) emission towards J1148 + 5251 
from: D-array only data imaged with a weighting parameter robust of 1 
(black), C&D-array data imaged with robust of 1 (blue), C&D-airay data 
imaged with robust of 0 (yellow) and C&D-an'ay data imaged with robust 
of -1 (pink). 


difficult because of the moderate signal-to-noise and non-ideal uv 
coverage, both due to a shorter-than-requested observation time on 
the full C-array. The combined C&D observations yield a synthes¬ 
ized beam with broad wings, depending on weighting. To further 
explore the issue of the CO (7 = 2 —> 1) emission size, we carried 
out two tests. 

First, we looked at point-source line profiles at various resolu¬ 
tions. We considered the position of the peak emission from the D- 
array only map (derived from a Gaussian fit) as the position of the 
point source and took line profiles from the D-array only and C&D- 
array data imaged with three different weightings (robust = 1,0 
and -1). These spectra are shown in Fig.|^ Moving from a close- 
to-‘naturaT (robust = 1) weighting towards a close-to-‘uniform’ 
(robust = -1) one there is a clear decrease in the peak emission, 
as it drops from ~ 0.30 to ^ 0.15 mJy beam^'. Fitting Gaussians 
to each of the spectra, we find a difference of ~ 60% between the 
D-array only spectral amplitude and the C&D-array robust = -1 
spectral amplitude. The restoring beams for each of these maps are 
listed in Table [T] Figure]^ shows that the peak flux decreases with 
increasing resolutions, and hence the source is resolved at higher 
resolution. 

Second, we imaged the combined C&D-array dataset. This is 
problematic as well due to the different weightings required for op¬ 
timal signal-to-noise by data from the two arrays and the broad line 
wings resulting from having significantly more time on-source with 
the D array versus the C array. To avoid these issues, we CLEANed 

the integrated line emission (velocities between-216 km s“’ and 

~ 189 km s“‘) to a depth where the synthesized beam’s broad wings 
are removed and explored different Briggs weights. Figure|^shows 
contours of the J1148+5251 velocity-integrated CO (7 = 2 —> 1) 
emission imaged with a robust of 0.25 overlaid on to the [C ii] emis¬ 
sion observed by jCicone et al.|j2015) . The spatial resolution of the 
CO (7 = 2 —» 1) observations is T.'07 x 0'.'97, just slightly higher 
than that of the [C ii] observation, as shown in the bottom left comer 
of the map. 

The main emission region is less than 1 arcsec in extent, al¬ 
though there is some diffuse emission to the south west of the 
source at 3cr which needs further confirmation. A formal elliptical 
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Figure 4. Velocity-integrated VLA (C & D configurations) contour map 
of the CO {J = 2 —> 1) line emission towards JI148 + 5251 overlaid 
on to [Cii] emission observed by |Cicone et al.H2015t . The [Cii] emis¬ 
sion is integrated over a larger velocity range that the CO one. The cross 
marks the position of the optical quasar jWhite et al.|2005y Contours are 
at (-2,-1.41,1.41,2,2.82,4,5.64,8)xo- (Icr = 25pJy beam-'). The syn¬ 
thesized beam of the CO data is I'.'OV X (y.'97 and that of the [C ii] data is 
1'.'33 X 1'.'23. 


Gaussian fit to the robust = 0.5 image returns a deconvolved major 
axis of O'.'89 ± O'.'46, and an undetermined minor axis and position 
angle. So we marginally resolve the emission in one direction and 
measure a size consistent with that found by |Riechers et al.H2009^ 
for the CO (7 = 7 —> 6) emission. 

The peak flux of 157 ± 17 pjy beam-' is measured close to 
the position of the optical quasar (see cross in Fig. [White et al.| 
|2005| > and agrees very well with the position of the [C ii] emission. 

Much deeper observations are needed to determine the source 
morphology as a function of velocity. 

We do not detect any continuum in the CO (7 = 2 —» 1) 
observations down to a limit of Icr = 6.7 pJy. 


4 ANALYSIS 


The line width of 298.05 ± 25.87 km s-' we derive for the CO 
(7 = 2 —» 1) emission agrees well with the width measured in 
previous studies Pert oldi e t al.|2003b{|Riechers et al.|2009 [ |Walter | 
|et al.|200^ . Hence, the low and high order transitions have similar 
line profiles in J1148 -l- 5251. 

The FIR-to-radio spectral energy distribution fit from|Wang et| 
[^ ( [2008) uses an emissivity index p = 1.6, a dust temperature T = 
55 K I Beelen^etalJ200^ and a FIR-to-radio ratio q = 2.34 ( |Yun,| 
[Reddy & Condon|2001| ) and predicts an observed flux density of ~ 
9 pJy, with a large uncertainty depending on q. Our non-detection is 
just below the predicted value. This result confirms a/3 value that is 
at least as steep as 1.6, and possibly steeper, thereby strengthening 
the argument for little cold dust in J1148 -l- 5251. 


4.1 Gas mass estimate 

Far infrared CO luminosity ratios and linewidths in quasars at high 
redshifts are similar to what is observed in local ultra-luminous in¬ 
frared galaxies (ULIRGs) jSolomon & Vanden Bout|2005[[^rilli[ 


[& Walter|2013[ l so we adopt the CO luminosity to H 2 mass con¬ 
version factor of or = 0.8 Mq (K km s-' pc^)-' used for ULIRGs 
( [Downes & Solomon[I998[ l. From [CarilIi & Walter| ( [20I3[ l, the line 
ratio Lpoz-i /^coi-o typically very close to unity for high redshift 
galaxies, so we can approximate L^oi-o 3.2 ± 0.3 x 10*° K km 
S-* pc^. Thus, from the CO (J= 2 —> 1) luminosity we derive a total 
molecular gas mass of Mh^ = 2.6 x 10*° Mq. 


4.2 Dynamical mass estimate 


Assuming a rotating disc geometry for the line-emitting gas in 
J1148 + 5251 we can estimate the dynamical mass of the emitting 
region using Mdy„ ~ 1.16 x lO^v^j^^, D Mq, where D is the dia¬ 
meter in kpc we measured for the emission and Vcirc is the FWHM 


in km s-* of the line profile. Following Wang et al. (2013i we 
use the approximation Vdrc sin i = 0.75 FWHMco(y= 2 ^i)- Tak¬ 
ing the size of the source as 8.1 kpc (the major axis) implies 
Mdyn sin^ i = 5.3 x 10*° Mq. For an inclination angle of 63° 
as found from the ratio of the minor to major axis, the dynam¬ 
ical mass is Mjyn = 6.7 x 10*° Mq, in line with previously de¬ 
rived values. Walter et al. (2004 > quote Mjyn = 4.5 x 10*° Mq (or 
Afdyn = 5.5x10*° Mq if accounting for an inclination of 65°) within 
a disc of diameter 5 kpc. [Riechers et al.[ j2009^ use the [C ii] line 
emission to find a dynamical mass of Mdyn = 1.5 x 10*° Mq for a 
region with radius 1.5 kpc. It has to be noted that there is substan¬ 
tial uncertainty to modelling assumptions (rotating disc), and model 
parameters (size, inclination) in our dynamic mass estimate. High 
resolution imaging in [C ii] or CO at much higher signal-to-noise is 
required to better understand the galaxy dynamics in 11148 - 1 - 5251. 


4.3 CO Line Excitation 

We have expanded the large velocity gradient (LVG) model of the 
line excitation in J1148 - 1 - 5251 calculated by [Riechers et al.[j2009[ l 
in which the kinetic gas temperature and density are considered 
as free parameters. A fixed 3 : 1 H 2 ortho-to-para ratio was used 
due to the relative statistical weights of the symmetrical and anti- 
symmetrical eigenstates of the wavefunction. This model uses the 
[Flower| ( |2001^ CO collision rates and a cosmic microwave back¬ 
ground temperature of 20.25 K at z = 6.42. We kept the parameters 


adopted by [Riechers et al.] ([2009[l: a CO abunda nce per velocity 


gradient of [CO]/(dv/dr) = 1 x IQ-^ pc/(km s-*) (WeiB, Walter & 

|Scoville|2005l|Riechers et aL|2006 

[WeiB et al.|2007ail and a CO 

disc radius of 2.5 kpc (jRiechers et a 

.|2009[ Walter et al.|2004[(. 


The CO (7 = 2 —> 1) data point fits well on the curve de¬ 
scribed by the best solution (Fig.|^. This was obtained for a spher¬ 
ical, single-component model with a CO disc filling factor of 0.16, 
Tun = 50 K, andpga4H2)= lO'* ^ cm-\ 


5 DISCUSSION 

We have observed the CO (7 = 2 —» 1) transition towards 11148 - 1 - 
5251 and find supporting evidence that 11148 - 1 - 5251 is similar to 
lower redshift quasars in terms of its high CO line excitation and 
compact size. We marginally resolve the main emission region to 
have a major axis of 0'.'89. We also detect tentative evidence for 
diffuse emission extending to the south west. 

Considering the line excitation in J1148 - 1 - 5251, our data are 
used to provide further constraints on previous LVG models, res¬ 
ulting in average physical conditions that are similar to those ob¬ 
served in other high-redshift quasars ([Weifi et al.|2007bl [Riechers] 
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Figure 5. CO excitation ladder (symbols) and LVG model (line) for J1148 + 
5251. The upper limit for CO (/ = 1 ^ 0) and the (continuum corrected) 
data point for CO (J= 6 —* 5) are taken from |Bertoldi et ar]j20Q3a). The 
CO (7 = 7 —> 6) (continuum corrected) data point is taken from |Riechers| 
|et al.|p009| an the CO (J= 3 —> 2) data point (corrected for missing flux 
in the line wings) is taken from [Walter et al^j20Q3| . The square marks 
the data point from the current observations. The model predicts kinetic 
temperatures and gas densities of Tjun = 50 K andpgas(H 2 )= lO"^-^ cm“^. 


|et al.|2009) . We find Ttin = 50 K extending from low to high order 
transitions. As such, it is possible that the emission in the low order 
transitions of CO may be associated with the same warm, highly 
excited star-forming gas traced by higher-i line transitions. This 
would support previous claims that the mid-T CO lines {J< 3) may 
also be good indicators of the total amount of cool molecular gas 
in high redshift quasars ( [Riechers et al.|2006l|2011^ . 

In the local Universe, the traditional view of ultraluminous in- 
fared galaxies (ULIRGs) and AGN has been that these objects rep¬ 
resent different phases in an evolutionary sequence ( [Sanders et al.| 
|1988| l. A heavily dust obscured galaxy undergoing vigorous star- 
formation will be observed as a ULIRG, while a buried AGN con¬ 
tinues to grow a supermassive blackhole. This may eventually give 
rise to a far-infrared luminous quasar, followed by an optically lu¬ 
minous quasar with only a modest gas reservoir remaining. This 
scenario may be similar to what is observed in high-redshift quas¬ 
ars, like J1148 -l- 5251, which may have gone through an earlier 
phase of starburst activity before the AGN had formed. Exten¬ 
ded reservoirs of cold molecular gas have been detected via CO 
line emission in z - 2 SMGs ( jCarilli et al.|2010[ ?; |Ivison et al.j 
[201 1[ [lechers et al.|2010||201 Ibjc) , suggesting that they may be 
in the process of forming through gas-rich mergers. A population 
of z ~ 6 SMGs has recently been discovered (e.g. Riechers et al. 
2013; Weiss et al. 2013), so it is possible that these may represent 
an early phase of the luminous z ~ 6 quasars such as J1148 -(- 5251. 
These objects are likely to evolve into the most massive early-type 
galaxies observed in the local Universe. However, in the case of the 
single known z ~ 6 SMG, the lack of a hot dust component detec¬ 
ted in the Herschel bands makes this object inconsistent with such 


a scenario. More work is required to establish the validity of the 
z ~6 SMGs and z ~ 6 quasars connection. 

We calculate new estimates of the gas and dynamical mass 
from the CO (7 = 2 —> 1) line luminosity. Our results agree with 
previous measurements, but are less prone to biases due to high 
order molecular line excitation. The gas mass we estimate is ~ 63% 
of the dynamical mass, suggesting a low value for a. This supports 
a departure from the Mbh - Cbuige relationship, in that the of 
the host galaxy is well below that expected from the low z relation. 
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